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^ . We study the gravitational collapse of an overdensity of nonrelativistic matter under the 

^ i action of gravity and a chameleon scalar field. We show that the spherical collapse model 

'^ I is modified by the presence of a chameleon field. In particular, we find that even though 

Q ' the chameleon effects can be potentially large at small scales, for a large enough initial 

size of the inhomogeneity the collapsing region possesses a thin shell that shields the 
modification of gravity induced by the chameleon field, recovering the standard gravity 
k> ! results. We analyse the behaviour of a collapsing shell in a cosmological setting in the 

H [ presence of a thin shell and find that, in contrast to the usual case, the critical density 

for collapse in principle depends on the initial comoving size of the inhomogeneity. 
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1 Introduction 

Fundamental scalar particles have not been found in Nature yet. Many examples have 
been postulated over the last forty years; the most illustrious being the Higgs boson, a 
cornerstone of the Standard Model (SM), which may well be the first one to be discovered 
at the Large Hadron Collider. The existence of scalar fields has also been invoked in 
order to explain the inflationary phase in the early universe [1]. In the same vein, the 
late time acceleration of the expansion of the universe could result from the presence of 
dark energy scalar fields[2, 3, 4, 5]. In all these cases, the coupling of scalar particles to 
normal matter leads to new interactions and possibly variations of parameters such as 
masses and coupling constants [6]. 

The physics of scalar fields coupled to matter depends crucially on the scalar field 
mass. For instance, the Higgs boson is expected to have a large mass (m^ > 115 GeV) 
compared to the Hubble scale today {Hq = 10"^^ GeV) implying that the new interaction 
mediated by the Higgs boson is of very short range. Moreover the Higgs field has been 
frozen at the minimum of its potential since the early stages of the universe. On the other 
hand, several extensions of the Standard Model predict the existence of light scalar fields. 
Some of these field might be moduli of string compactifications or the dilaton. Another 
class comprises the fields which can appear as a consequence of the spontaneous symmetry 
breaking of a global quasi-exact symmetry. These pseudo-Goldstone bosons may play a 
relevant role in explaining the recent stage of acceleration observed in cosmological data. 
In this case, the scalar fields are nearly massless as their mass should be of the order 
of the Hubble rate now. All these light scalar fields can couple to matter and hence 
introduce new long range forces which would appear as a modification of gravity. The 
new interaction would also induce a non-conservation of the matter field stress-energy 
tensor, which can be interpreted as arising from mass-varying matter fields. Limits on 
this type of interactions arise from SNIa and CMB data [7]. As we will see below, 
alternatively one can choose to work in the so-called Jordan frame, where the matter 
field stress-energy tensor is conserved. In this case, Newton's constant varies and its 
variation since Big Bang Nucleosynthesis should not exceed some forty percent [6]. 

Modifications of gravity should not be taken lightly. There are stringent limits on the 
couplings of scalar fields to normal baryonic matter arising from the search of a fifth force 
or deviations from the usual gravitational force at different scales, from the laboratory to 
the solar system. In particular, the Cassini bound [8] implies that the coupling to matter 
should be less than 10~^. Several mechanisms have been proposed to alleviate this prob- 
lem. The first one which applies to the dilaton in the strong string coupling limit is due 
to Damour and Polyakov [9] and it asserts that the field is attracted by the zero coupling 
limit cosmologically. Another mechanism occurs in the Dvali-Gabadadze-Porrati model 
[10] whereby the brane bending mode, whose presence would lead to strong deviations 
from Newtonian gravity, is suppressed at short distance thanks to the Vainshtein mecha- 
nism. Another possibility is to generate a mass term for the scalar field which depends on 
the matter density of the environment. This is implemented in the so-called chameleon 
mechanism [11, 12], where a coupling between a scalar field (the chameleon field) and 
matter arises from a conformal rescaling of the metric in the matter sector. The nonlinear 
interaction of the chameleon field generates a shielding, known as the thin-shell effect. 



that diminishes the modifications of gravity for dense and large enough bodies. 

In this paper we are interested in modifications of General Relativity in chameleon 
models, focusing on the possible consequences for the growth and collapse of large scale 
structures in the universe. The growth of perturbations in chameleon models in the linear 
regime was studied in [13] . Recently, there has been a new bout of interest in studying 
the infiuence of modified gravity models on the formation of large scale structures. Indeed 
most of the aforementioned models of modified gravity show very little deviation from 
a pure ACDM model since Big Bang Nucleosynthesis. Differences only emerge at the 
perturbative level in the way structures grow. In that respect it is of high interest to 
go beyond the linear level and study the nonlinear regime of structure formation. A 
semi-analytical tool that has been used to analyse the nonlinear stages of gravitational 
collapse is the so-called spherical collapse model [14]. 

The spherical collapse model has been recently used in models with a simple Yukawa- 
type modification of gravity [15], in the so-called f{R) models [16], in braneworld cos- 
mologies [17] and in models which allow for dark energy fiuctuations [18]. We will focus 
on the collapse of objects with a thin shell that arises in the context of chameleon models. 

This paper is organized as follows. In section 2 we give a brief introduction to 
chameleon models and discuss the modifications of gravity which occur in these models, 
both for geodesies and the growth of perturbations. In section 3 we solve the dynamical 
equation for the scalar field in a spherical configuration of matter overdensity, finding its 
profile and studying the effects of the existence of a thin shell which effectively shields the 
modifications of gravity. These profiles serve as templates for the chameleon behaviour 
during the spherical collapse. Section 4 presents an analytical discussion of the conse- 
quences of the chameleon models for the spherical collapse models of large scale structure 
formation. Phenomenological consequences are expounded in section 5 with a simple 
example before presenting our conclusions. 

2 Modification of Gravity with Chameleons 

2.1 Chameleon models 

The action governing the dynamics of the chameleon field is given by 

S = J d'xV^l^n-^g^''d,(l>d,(l>-Vi(l>)\+ J d'xy^CUi^l:},~9,u), (2.1) 

where Qfj^i, is the metric in the Einstein frame, g^^ is the metric in the Jordan frame, 
Mpi = (StiGn)^^^'^ = k'^^ is the reduced Planck mass, TZ is the Ricci scalar for the 
Einstein-frame metric, and ipm are various matter fields labeled by i. The Einstein-frame 
metric g^j^iy and Jordan-frame metric g^^ are related by the conformal rescaling 

~g,u = e'^^/^-^g,u = n\<t>)g,,, (2.2) 

where /3 is a dimensionless constant [19]. In chameleon models this constant is assumed 
to be the same for all types of matter, respecting the weak equivalence principle where 
all material particles follow the same metric. This conformal coupling of with matter 
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particles is a key ingredient of the model, and as a result the excitations of each matter 
field ^jJm follow the geodesies of the Jordan frame metric. As a consequence, the physical 
energy momentum tensor of matter T)i™ is defined in the Jordan frame [21, 22], where 
it is covariantly conserved 

y^^{m)M^ = 0. (2.3) 

It is related to the energy- momentum tensor in the Einstein frame, T^™ , by T ^^ = 

In this paper we restrict our attention to non-relativistic matter for which the only 
non-vanishing component of T ^y is T q = — p, giving an energy density of matter 
p = Q'^p in the Einstein frame. In the Einstein frame, the equations of motion following 
from Eq. (2.1) are then 



^ ^u 



Ki 



T^T^ + d,<PdA - ld''<pd^<pg,u - Vg^, 



(2.4) 



where T*-™-* = g^^T^u^ = — pfi^. Now consider a spatially fiat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) geometry with Einstein frame metric 

ds^ = -dt^ + a{tfdxidx\ (2.6) 

where t is cosmic time. Then, for nonrelativistic matter and a spatially uniform scalar 
field, Eqs. (2.4), (2.5) and (2.3) reduce to 



PI 
a 1 



a 6M2j 



pp + p^ + 3p<^] (2.^ 



+ 3/70 = — ^-fip (2.9) 

dcp 

'p+3Hp = 0. (2.10) 

Here Q' = dQ/d(f), a dot denotes a derivative with respect to cosmic time t, and the energy 
density of the chameleon field is p^ = 0^/2 + V whilst its pressure is p^ = 0^/2 — V. The 
quantity p is defined by 

p = n-^p (2.11) 

and is conserved as a consequence of (2.3) (for non-relativistic matter only). Equivalently 
the Einstein frame energy density p satisfies 

p + 3Hp = (3k4P^. (2.12) 

It is important to draw a distinction between the three energy densities p, p and p. The 
latter is the Jordan frame energy momentum which is always conserved, independently of 
the equation of state of matter. As a result, the density p will appear in the equations for 
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conservation of matter and for the fluid velocity. The Einstein density p is the source of 
the Newton potential as it appears in the Poisson equation. For non-relativistic matter, 
the chameleon field dynamics are influenced by the conserved matter density p. It appears 
as a source term which modifies the effective potential of the chameleon. 

Indeed from the Klein-Gordon equation one can define an effective potential for the 
chameleon,"^ 

VM) = \/(0) + pn = \/(0) + pe^'t>l^^\ (2.13) 

It follows that V^ff has a minimum at (pmin{p) which is solution to 

'"' ' -p. (2.14) 



M, 



PI 



Notice that the source term is the Einstein energy density which depends on the chameleon 
field via the conformal factor Q, since p = Qp. The mass at the minimum is given by 






^ =-77¥ + T?rP- (2-15) 



PI 
As ^ < and ^^ > for successful chameleon models [11], 

m^ > 3I3^H^ (2.16) 

where we have assumed that the chameleon energy density is negligible compared to p, 

P<t> <P- 

In contrast to usual quintessence models, in most cases the mass of the chameleon is 
much larger than the Hubble parameter. An explicit example will be given in section 5. 
Thus it is possible to have fluctuations in the chameleon field at scales much smaller than 
the horizon scale which can be of cosmological importance for structure formation. 

The minimum of the effective potential is an attractor, and the chameleon therefore 
settles down there, and evolves slowly since at least Big Bang Nucleosynthesis. In that 
case, the potential energy of the chameleon is nearly constant and the model is equivalent 
to a ACDM model. We will see that the dynamics of chameleon models are not equivalent 
to a ACDM model at the level of perturbations and when the overdensities experience a 
spherical collapse phase. 

2.2 Geodesies 

We are interested in the dynamics of CDM in the presence of a chameleon field back- 
ground. The trajectories of the CDM particles satisfy the geodesic equation in the Jordan 

frame 

cPx'^ ~ dx" dxP 

where 

df"^ = -g^^dx^'dx'' (2.18) 



*Here we assume that, since p is conserved, dp/dcf) — 0. 



and the Christoffel symbol are calculated with the metric g^y. In the following we consider 
the perturbed Jordan frame metric (in conformal-Newtonian gauge) 

^§2 ^ ^21^^-^ ^_i^-^ ^ 2<j)N)drf + (1 - 2<j)N)dx'dx,] , (2.19) 

where t] is the conformal time and the Jordan frame scale factor a is related to the 
Einstein frame scale factor a by 

a = rta. (2.20) 

The scalar metric perturbation (pN is related to the perturbation in the Newtonian grav- 
itational potential and is sourced by the CDM density perturbation S defined through 

p = Poo{l + 5) (2.21) 

where poo is the background energy density. Notice that when <^ Mpi, the Jordan 
frame metric can be expanded as 

dS^ = a^{r]) [-{1 + 207V + 2(5K^(t))drf + (1 - 20^ + 2/3/t40)rfa;Va;i] . (2.22) 

Hence, in the Jordan frame in which particles couple to gravity, two distinct Newtonian 
potentials appear whose difference comes from the chameleonic field. We will see that 
the existence of a chameleonic contribution to Newton's potential identified from the drf 
element affects the motion of particles along geodesies. 

Before continuing, we now specify certain approximations which will be made through- 
out the remainder of this paper. First, we assume that 

07V < 1 (2.23) 

during structure formation. Then we assume that: 



u 



<P'n « n, j^-^-H (2.24) 

iWpi 

where ' = d/drj and "H = dlnd/drj. This implies, for example, that the Hubble parameters 
in the Einstein and Jordan frame are essentially identical: "H = 'H + /30'/Mpi ^ "H. It also 
implies from (2.12) that p, as well as p, are approximately conserved. This is compatible 
with our next approximation, namely that 

< 1 =^ n{(j)) ~ 1. (2.25) 



Mpi 

Finally, in the non-relativistic limit, all spatial gradients dominate over time variations, 
so that for example 

|<9^0| > 0'. (2.26) 

Using (2.23) and (2.24), the geodesic equation in the Jordan frame, Eq. (2.17), can 
be simplified. From (2.22) the Christoffel symbols are given by 

Fqq = 529*^, fjo ^ nS] (2.27) 
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where 

^ = 0^ + /3K4(f). (2.28) 

Hence the geodesic equation becomes 

rfV ndx' ^. , , 

-^^ + 2--^ = -d'^. 2.29 

dr"^ a dr 

where we have used that 

dx^ = dr] = -df (2.30) 

from the zero component of the geodesic equation. Thus in proper time and in comoving 
coordinates, particles follow Newtonian trajectories under the influence of the effective 
gravitational potential \1' given in Eq. (2.28). In other words, matter moves under the 
influence of a new Newtonian potential modified by the profile of the scalar field. 

2.3 Perturbations 

Chameleon theories are very different from a cosmological constant at the perturbative 
level. This can have important effects on structure formation, as we will see shortly. 

The Lagrangian description of trajectories is equivalent to an Eulerian one in which 
the velocity field f * = ^ obeys the Euler equation. Let us define the unit 4-velocity 



drj 

dx'^ 



U 



fJ- 



{u^u^~g^u = -1) (2.31) 



v" - Uv' + v^djv' 



(2.32) 



dr 
so that in the non-relativistic limit, -u* = v^ /a. Then using the chain rule, we obtain 

d'^x' _ 1 
df^ a^ 
Thus the geodesic equation reduces to the Euler equation 

v" + v^djv' = -Uv' - 5'^dj^ (2.33) 

for the nonrelativistic fluid. Notice the presence of the effective Newton potential ^ . 
The divergence of (2.33) gives 

e' + ?/0 + — = -A^ (2.34) 

o 

where = (9jf*. Furthermore, energy conservation in the Jordan frame implies that the 
density perturbation 5, defined in Eq. (2.21) satisfies 

(5' + 0(1 + (5) = 0. (2.35) 

Combining eq.(2.34) and eq.(2.35) leads to a nonlinear equation for the evolution of 
perturbations: 

<5" + m'-|-^ = (l + 5)A^. (2.36) 

3 1 + 

We will use this equation in order to study the spherical collapse. 
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2.4 Modified gravity 

For the following discussion it will be more transparent to use physical coordinates 

r' = ax' and V^ = d/dr\ (2.37) 

Then the geodesic equation Eq. (2.29) is rewritten as 

(fr^ . ^ l(Pa j 

dr'^ a dr"^ 

. / 1 (i^a 2 

= -V ^ r 

V 2Srff2 

= -V^^^. (2.38) 

This is just Newton's law in an expanding universe, and matter moves under the influence 
of the modified Newtonian potential '^ n. On using Eq.(2.28), 

'^N = ((l)N-—:T^rA+f3K4(p 
\ 2a dr^ J 

= $jv + /3/t40- (2.39) 

The 'total Newtonian' potential $Ar satisfies the Poisson equation sourced by the total 
density of matter. Indeed, let us now work under the assumptions given in Eqs. (2.23)- 
(2.25). Then from the zero-zero component of Einstein's equations (2.4) in the Einstein 
frame, 0Ar is solution to 

A0^ = ^Poo<5, (2.40) 



whereas we also have 



1 d'^a a AnGr 



(Poo + (l + 3w^)p<^). (2.41) 



a dr^ a 3 

Thus from Eq. (2.39) it follows the Poisson equation reads 

2 2 

A<l>^ = ^ [p + (1 + Sw^)p^] ^ ^p (2.42) 

where we can neglect the chameleon energy density until very recently in the history of 
universe. 

The resulting modifications of gravity are drastically illustrated for models with V = 
in Minkowski space and p = mQ5^^\ namely a massive point particle. In this case the 
Klein-Gordon equation (2.5) reduces to 

A<p = /3n,moS^'^ =^ = _^^ (2.43) 



Thus from Eq. (2.28) 



Attv 



^ = - (1 + 2/3^) —mo. (2.44) 



p 



oo 




^b R 



Pc ^'^'^ 0-W 



-*r 



{ 



Figure 1: Scalar field profile with no thin-shell in the presence of a spherical overdensity 
with constant density pc and radius R immersed in a background of density poo- 

Hence there is a modification of Newton's law depending on the coupling /3 

GN,eS = 7Gn, (2.45) 

where 

7 - 1 = 2/32 (2.46) 

Of course, large /3's would be excluded by the Cassini bound on the existence of 
fifth forces. This is not the case for chameleon theories. Indeed two mechanisms are at 
play. The first is that the chameleon develops a density dependent minimum, and at 
the minimum its mass increases with the density of the environment. This is sufficient 
in the atmosphere where the density is large, and this implies that there is no effect on 
Galileo's Pisa tower experiment. In a sparse environment such as the solar system or a 
laboratory vacuum chamber, a second mechanism must be invoked: the existence of a 
thin shell which implies that 7 — 1 ^ 2/3^ for massive bodies such as the moon, the sun 
and laboratory test bodies. 

3 Chameleon Profiles in Spherical Matter Configu- 
rations 

We now find approximate solutions of the Klein-Gordon equation Eq. (2.5) in physical 
coordinates. More precisely, in this section we assume that the matter source is a perfect 
sphere of constant radius R with constant homogeneous density p^. It is immersed in a 
homogeneous cosmo logical background of density poo *^ Pc(see figure 1). On sub-horizon 
scales where one can neglect time derivatives in the evolution equation of the chameleon, 
one can treat this case as a quasi-static problem. Our analysis is a generalisation of that 
presented in [11]. 

For the following discussion it will be useful to define the field values 0c and 0oo 
as those which extremise the chameleon effective potential for the densities pc and p^c, 
respectively. Hence from Eq. (2.14), 

Koo = -7- = -/3«4Pc,oo- (3.1) 



3.1 Small bodies, no shell 

First we assume that the sphere is small (in a sense which will be made precise below) 
and its Newton potential at the surface is small also. Let (j)b denote the value of the field 
at the core. We expect (f)b to differ significantly from (pc and indeed that it should be 
closer to 0oo in the cosmological background. 

Inside the overdensity, the Klein-Gordon equation Eq. (2.5) reduces to 

rf20- 2d(t)~ dV ^ , , 

ar-^ r ar dcp 

We expect 

0"(r) = 0b + 5(t)-{r) with 50" (r) < 0^ (3.3) 

where the value of the field at the core, 0f, = 0^(0), must be determined. Substitution 
into (3.2) and expanding to first order gives 

%.l'^^<(r-,.Hv^-K (3.4) 

where Vl = d(f,V((j)b) and ml = V"{(f)b)- Thus the solution inside is simply 

r^4,„ + Yl^(i-?^^:hi!iiL) (3.5) 

m^ \ rribr J 

where the two integration constants have been fixed through the conditions d(f)' /dr = 
at r = and 0~(O) = 0^. 

Outside the overdensity, 0+ = 0oo + Sep, and following the same procedure leads to 
the equation of motion 

d'<P+ 2d<P+ , 

-7^ + --!— ="^oo(0 -0oo) 3.6 

dr"^ r dr 

where m^ = V"{(f)oo), with solution 

r)f,-moo{r-R) 

0+ = 0OO + — . (3.7) 

r 

Here D is another integration constant and we have imposed that 0+(r — ?> oo) = 0oo- 
The two solutions must match on the surface of the overdensity, namely 



0+(i?) = 0-(i?) 



dr 
This leads to 



R dr 



(3.8) 



R 



and 



■'oo 



VI - v 

D = — T-, ^-——(sinhmbR — rribR cosh m{,R) (3.9) 

m^(l +m^R) 



Vl — V^f sinh m^i? 1 sinh m^-R — m^-R cosh mfei? i /'q m^ 



ml \ wibR (1 + moo-R) mibR 
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In most interesting applications, the overdensity is small compared to the Compton wave- 
length m^. 

On assuming that both m^R <^ 1 and rrihR ^ 1, these give 

D = ^(K'-K) (3.11) 

0. = 0oo + ^(V;'-n')^'- (3.12) 



This second equation determines (pb- on searching for a solution of the form 0^ = 0oo+c 
where 6(f)oo <^ 0oo we find 

(3 .._^o2 



■'oo 



iAp)R' (Ap = p,-poo). (3.13) 



"" "^ 3Mpi 
Using the condition 0^ < 0b we find that this description is valid provided 

«:4|0c-0oo|>2/3|0jv(i?)| (3.14) 

where 0Ar(-R) = —Gn^M/R is Newton's potential at the surface of the overdensity. We 
will call this the no-shell condition. 

The constant D is determined from the boundary conditions to be D = /3/t4(AM)/(47r) 
where AM = 47r(Ap)i?'^/3. Thus from Eq. (3.7), the potential outside reads 

^40"'' = fi;40oo — 2/3 . (3.15) 

r 

As a result, from Eq. (2.28) the contribution of the chameleon to the force acting on tests 
particles is given by 

^ = /3«:40oo - (1 + 2/3^)^^^^ (3.16) 

r 

The first term has no gradient, the second term gives the force on tests particles. It only 

involves the overdensity and its mass AM. This corresponds to an increase of Newton's 

force by (1 + 2/3^) 

7-1 = 2/3=^ (3.17) 

as obtained earlier for a test mass. 

When discussing the gravitational collapse of an initial overdensity, if the density 
perturbation is small and the size of the overdense region is not too large, then the inside 
of the collapsing sphere will be described by the results in this section. This is true as 
long as (pb does not converge towards 0c, i.e. when the no-shell condition is satisfied. 
When the no-shell condition is violated, the behaviour of the chameleon profile changes 
drastically. 

3.2 Large bodies, thin shells 

For large bodies, we expect the value of the field at the centre of the spherical body to 
be 0c, i.e. when Newton's potential is large enough, the solution is identically constant 
inside a shell 

(l)-{r) = (l), (0<r<i?,) (3.18) 
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(see Eq. (3.5) for 0^ = 0c)- Following [11], in the shell 

(r) = 0c + ^— I y + — 1 ^ {Rs<r<R) (3.19) 

while outside the body 

where we have not assumed anything about the magnitude of rriooR- Notice that the field 
outside the body is only determined by the mass inside the shell. The radius of the shell 
itself, Rs is determined from continuity of the fields at R, namely 0''(-R) = 0"'"(-R), or 

4 - ^» = r.Pp.%1 - Y^) + l-.PP.RHi + ^^) - I^^PpX (3.21) 

When rriooR ^ 1, this reduces to 

0,-0^ = ^(i?,2-i?2). (3.22) 

In the following we will always assume that rriooR <^ 1 which will correspond to the 
astrophysical situation of interest. 
Let 

Rs = R-AR (3.23) 

Then to first order in AR/R, Eq. (3.22) reduces to 

A^ = J^_Z^ (3.24) 

where $7v(-R) = —GnM/R is the total Newtonian potential at the surface of the body 
(see Eq. (2.42)) and M = AiiR^pc/S is the mass of the body. Eventually when the radius 
of the compact body is large enough, the shell appears. More precisely, this occurs at a 
radius Re which, from Eq. (3.24), is determined by 

«:4(0oo-0c) = 6/3|<l',v(i?c)|. (3.25) 

When (0OO - 0c) < 6/3Mpi$Ar (i?c), the shell is really thin. 

The modification of Newton's law outside the compact body is given by 

GnAM ^^2/ R^AGnM ,^^^, 

"^ = Kil3(l)^ - -^^ 2/3Ml--^j— ^. (3.26) 

Notice that the modification of gravity is entirely due to the mass inside the shell. The 
rest of the mass of the spherical body is completely shielded. The Laplacian of Newton's 
potential \1' at the surface of the overdensity is 

A^l^^^ = 47rG'jv(Pc - Poo) + StiGnP^ (l - ^] p, (3.27) 



\r=R -■•"it\rc roo/ i ivr- i - ^g 
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whilst the total Newtonian potential "^n including the contribution of the cosmological 
density inside the sphere, see Eq. (2.38), is given by 

A^ifNl=R = 4vrG;vPc + SttGjv/?' U - ^^ Pc- (3.28) 

This is the Newtonian potential multiplied by 

7(i?,i?s) = l + 2/32(^l-||^ (3.29) 

extrapolating between General Relativity {'-/ = 1, Rs = R) and a small body regime 
(7 = 1 + 2/32, ji^ = oy 

4 Spherical Collapse 

In the previous section we have found a description of the modification of gravity for 
chameleon theories in case of a static, spherically symmetric, overdensity. We now gen- 
eralise this discussion to the dynamical situation in which the density inside and outside 
the body evolve on cosmological time scales. 

4.1 Initial conditions 

Before writing down the set of nonlinear equations governing the evolution of the system, 
it is very useful to analyse the astrophysical situation in which an initial inhomogene- 
ity appears as a spherical body of radius R immersed in a homogeneous background. 
Depending on the radius, the initial body may or may not have a chameleonic shell. 

The radius Re below which no shell is present is given implicitly in Eq. (3.25). This 
simplifies when pc = {1 + Si)poo with 5j <C 1. First of all the minimum of the effective 
potential is such that 

«:4(0c - 0oo) = -^^<5, (4.1) 



"^So 



and therefore from Eq. (3.25) 

fl. = ^. (4.2) 

rrioo 

For large objects with R > R^, we will use the thin shell solution^. In particular, this 
will imply that the gravitational collapse of these objects is not too strongly modified 
compared to the General Relativity case. 

4.2 Collapse 

Hence let us consider an initial spherical overdensity with a thin-shell, thus the initial 
radius Ri > Re so that initially the thin shell is between Rg < t < Ri. We now try to 
understand qualitatively the evolution of the system. As we shall see, the dynamics of 



^We always check that the condition rUooR <C 1 is verified. 
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the nonrelativistic particles inside the overdensity depend crucially on whether they are 
within the thin shell or not. 

Let g* denote the initial position of the particle. As shown in the previous section, 
see Eq. (3.18), if |g*| < Rs, the chameleon field is constant and the particle experiences 
no extra force compared to gravity 

7(r) = 1, r < Rs. (4.3) 

So initially, the particles on a spherical shell of radius |g*| < Rs collapse as in General 
Relativity. For particles on shells within the thin-shell region itself, the initial force they 
feel is enhanced compared to gravity and 

7(r) = l + 2/32('l-^V r>Rs. (4.4) 

This implies that the spherical collapse of these spherical shells |g*| > Rs is initially 
accelerated compared to the shells with |g*| < Rg. As a result, due to the relative 
acceleration felt by particles at different radii, the density inside the thin-shell will not 
remain uniform: at a given instant t and following the particles inside the shell in a 
Lagrangian way, the thin-shell density becomes 

p(r,t)=p(g)det(^). (4.5) 

Of course this density depends on the dynamics of the nonrelativistic particles and there- 
fore on the chameleon field in a bootstrap fashion. 

The radius of the thin-shell Rs, defined as the boundary where no scalar force is felt, 
starts moving too, thus Rs = Rs(t). For radii smaller than Rs(t) we have 

= 0c, r<Rs{t). (4.6) 

Mass conservation implies that the mass initially outside the thin shell with r < Rs{t) is 
conserved 

Ms = ^^Rlit). (4.7) 

with a dynamical evolution as in General Relativity. Inside the thin shell the chameleon 
field satisfies a Poisson equation as the potential contribution is negligible [11] 

A0 = /3/€4p(r, t), Rs{t) <r < R{t) (4.8) 

where both the inner radius Rs{t) and outer radius R{t) of the shell are time-dependent. 
The density inside the shell is not uniform now. Let us define the mass inside the shell 

M{r,t)=An dr'r'^p{r',t) (4.9) 



Gauss' law implies that 



[r,t) = (5K,j M^dr' + <l>, (4.10) 
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This allows one to evaluate the force due to chameleon on a shell of radius r and therefore 
determine the motion of the nonrelativistic particles and the density of matter in the thin- 
shell. 

Outside the body, Gauss' law allows one to write the outer solution as 

^^,,^.p.,Mmi ,4,u) 

Attt 

as long as rriaor ^ 1. Matching at the surface r = R[t) the inner and outer solutions 
gives immediately 

,_,,.,., r MM,,-,, .,^w). ,4.12) 

jRsit) 47rr'2 47iR{t) 

This equation relates the radius of the outer shell R{t) to the radius of the inner shell 
Rs(t) and p{r,t). Notice that these equations are the generalisation of the ones obtained 
in section 3.2 when the density inside the spherical body is both constant in time and 
spatially uniform. Here this is not the case anymore due to the collapse of the spherical 
body. 

The previous description of the spherical collapse with the development of an inho- 
mogeneous thin-shell due to the different accelerations felt by the nonrelativistic particles 
at different radii inside the thin shell is very difficult to solve numerically. It involves a 
boostrap approach as the chameleon field in the thin shell is determined by the matter 
density, which in turn is determined by the motion of the matter particles under the 
infiuence of the chameleonic field. In the following, we will not attempt to numerically 
solve this tough problem. We will work with the lowest approximation assuming that the 
density inside the shell remains uniform spatially. More precisely we will follow the same 
approximation scheme as in the gravitational collapse with no chameleon. We assume 
that there is no shell-crossing and therefore the spherical shells at ordered radii ri < r2 
remain ordered. This implies that the mass in each shell is constant in time. In this case, 
starting from a homogeneous distribution of matter we have 

M{r,t)=M'-^^ (4.13) 

where Ri is the size of the initial overdensity and M its total mass. As a result we find 
that 

M fr' , Rlit) 3^2, 



0(^, t) = ^^^^^ (,y + r ~ 2^'^^V ^ "^^ ^^'^^^ 

This corresponds to the solution with a uniform density pc in the shell as already obtained 
in section 3.2. This first order approximation will give us a reliable description of the 
collapse when the size of the thin shell is very small and inhomogeneities have little time 
to develop. Of course, a more thorough study of the full set of equations is required to 
assess the full validity of this approximation. In general, shells will cross in a finite time 
leading to the formation of caustics. This is left for future work. 

Moreover, when applying the results on the spherical collapse with a thin shell, one 
should be wary of the fact that in realistic cases for astrophysical overdensities, the pres- 
ence of inhomogeneities breaking an exact spherical symmetry will lead to even greater 
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difficulties. Indeed, the modification of gravity induced by the chameleon will be non- 
spherical implying that matter particles will feel a nonspherical chameleon force which 
could therefore amplify the initial inhomogeneity. A full numerical simulation of the 
dynamics of a thin-shelled and a slightly inhomogeneous overdensity is certainly worth 
carrying out. 

4.3 Sphere dynamics with a thin shell 

Assuming that the density inside the thin-shell remains uniform, we can write down a 
simplified equation for the dynamics of the spherical collapse. Using the Poisson equation 
modified by the interaction of matter with the chameleon field, Eq. (2.38), and mass 
conservation one gets 

f = -^ [P0(1 + ^w^) + liR, Rs)Pc] (4.15) 

with' = d/dt and 7 is given in Eq. (3.29). This equation includes the chameleon enhance- 
ment of gravity taking into account the thin-shell effect, where Rg is the radius of the 
shell. The normal gravity case is obtained either for /3 = (no coupling) or if the shell is 
very thin, Rs/R — )> 1. 

We will assume that pip is constant {wip = — 1) and write 

^^ = 4^ = R^ = ^i ' ('•'') 

where again we have used mass conservation, and the usual assumption Q ^ 1: 

4:71 

M = M, = —R^p,{l + 6,), (4.17) 

Pi = p^ = poa^^ is the initial background density. We define R(t) as the radius of the 
collapsing spherical region in units of the initial comoving radius Ri/af. 

m = 1^ (4.18) 

with R{ti) = ai. (Notice that, as explained above, we have assumed f] ?« 1, so that we 
can identify the energy densities in the Jordan and Einstein frames, and a with a.) 
Using 

pi? = ^ (4.19) 

we can write Eq.(4.15) in terms of R as: 



R__H^ 
R~~^ 



I {-2nf + t^{R,R,)n^Z\l + 5,)^ (4.20) 



where we will see below that the only explicit dependence on Ri arises in the 7 factor. 
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4.4 The linear regime 

We have seen that the chameleon dynamics and the existence of a thin shell modifies 
the collapse of an initial spherical overdensity. It is particularly interesting to analyse 
the evolution of such a configuration in the linear regime extrapolated to values of the 
density contrast where the linear approximation is not valid anymore. Indeed, within 
this approximation, we can also compute the value of the critical density, Sc{z), used in 
the Press-Schechter formalism for structure formation [23]. It is defined as the linearly 
extrapolated density contrast from an initial value such that collapse occurs at a given 
redshift z. Practically, one must first find the initial density contrast required for collapse 
at z = 0. Then this initial density contrast evolves using the linearized form of eq.(2.36). 
Let us first consider the linear evolution of an overdensity with a thin shell. The 
Laplacian of Newton's potential \E' at the surface of the overdensity is given in Eq. (3.27): 



A^\r=Rit) = ^7iG^p^6 + SttG^v/S' ( 1 - 7^ ) Poo(l + S) (4.21) 



-R; 

i?3 



In physical time t and from Eq. (2.36), the linear equation for 6 valid only in the thin 
shell regime is: 



5 + 2H5 - AttG 



N 



1 + 4/3^1-1 



PooS = AttG 



N 



2/3M1-I 



Poo . (4.22) 



For P = (no coupling) or Rg = R (complete shielding) one recovers the standard linear 
evolution equation. Let us analyse briefiy the structure of this equation. First of all 
there is a forcing term which arises since the chameleon couples to all matter inside the 
sphere, whereas only the perturbed Newtonian potential contributes to the growth of 
perturbations. Then the term linear in 6 receives two contributions in 2/3^ when Rs — ?■ 0. 
This follows from the 2/3^(1 + 6)'^ term in the full nonlinear equation. This equation is 
valid as long as the thin shell condition is valid, i.e. as long as Rg > 0. In particular, the 
thin shell is very thin when K4,\(pc — (poo\ ^ 6/3|<l>Ar(i?(t))| where <l>Ar(i?(t)) is total Newton 
potential of the sphere. When the sphere is small enough and/or sparse enough, the thin 
shell regime is no longer valid. In this case, the field configuration for the chameleon is 
different. In particular we have seen that as soon as K,4\(f>c — 0oo| > 2/3|0Ar(_R(t))|, where 
0Ar(-R(t)) is the Newton potential due to the overdensity only, the value of the chameleon 
at the centre of the sphere is not 0c anymore. In this regime, the Laplacian of \E' is simply: 

A^|,=ij(i) = AtcGnPooS + MGnP^Poo5 (4.23) 

where the first term comes from 0jv and the second from the chameleon. In this case, 
the linear equation reduces to: 

5 + 2H5 - A-kGn [1 + 2/3^] p^5 = (4.24) 

Notice that the forcing term has disappeared and the effective Newton constant is simply 
(1 + 2/3^) larger than Gn- This is nothing but the usual linear equation deep inside the 
Compton wavelength of the chameleon field. This is justified here as this equation is valid 
for radii R{t) <^ v5/moo- 
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There is a sharp difference between the two regimes. When no shell is present, the 
linear evolution equation does not depend on R and therefore the critical density contrast 
does not depend on the initial value i?j. This is not the case anymore in the thin shell 
regime. Indeed, the evolution equation is sensitive to the ratio Rg/R which depends on 
the initial size Ri. As a result, the critical density contrast becomes scale dependent. 
This type of behaviour leads to a moving barrier problem for structure formation. The 
study of the connection between moving barriers and chameleons is left for future work. 

5 Phenomenology 

5.1 Inverse power law models 

In the following we will focus on models defined by an inverse power law potential and a 
constant coupling to matter [11, 12]. These models are the original chameleon ones and 
are defined by the potential 

A4+n 

V{<P) = At + ^^ + ... (5.1) 

where we neglect higher inverse powers of the chameleon field. The effective potential 
has a minimum at 

and the mass of the field at the minimum is given in Eq. (2.15). For this model (assuming 
< Mpi), 

In a cosmological setting, the Hubble parameter H ~ y^p/Mpi so that m 3> \f^B.. Hence, 
the field quickly relaxes to the minimum of the potential and sits there for most of the 
cosmological history, at least since matter equality. The fact that the scalar field can have 
a mass much larger than H is one of the main differences between chameleon models and 
usual quintessence models. At the minimum of the potential, 

^e//(0mm) ~A^ + p. (5.4) 

This implies that the model behaves like a pure cosmological constant in the recent 
past. The matter contribution is simply the usual energy density entering the Friedmann 
equation. Hence the main difference between these chameleon models and a ACDM model 
will be at the level of the growth of structures only. 

Alternatively, one can also characterise the chameleon potentials by the evolution 
of the chameleon mass moo (a), and the coupling /3 with the Cold Dark Matter sector. 
In the appendix we show that a power law evolution of the mass moo (a) = vdiiajaiY 
with n = — (6 + 2r)/(3 + 2r) leads to a power-law potential such as Eq. (5.1) and we 
explicitly compute A as a function of vrii and r for this case. In the following we will utilise 
the inverse power law as a useful parameterisation of the chameleon potential valid on 
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cosmological scales starting from matter equality, when the matter density evolves from 
p ^ lO^-'^'^g/cm^ to p ^ 10~^^g/cm^. 

In order for the chameleon effects to be significant, the Compton length scale m^ in 
the cosmological background should be larger than the scale of the perturbation R. As 
we will study perturbations at the Mpc scale, we require the initial mass rrii at matter 
equality to be m^^ = (9(0.02Mpc). Since we typically have small initial perturbations of 
the order of 6i = 10^^ at Zi = 10^, we estimate Re = (9(10^^Mpc) from Eq. (4.2) and 
we are justified to use the shell solutions for scales larger than Re. In terms of comoving 
scales, we consider scales ranging from Ri/ai = 1 to 100 Mpc for which the condition 
TJiooR <^ 1 is satisfied. Moreover, a thin shell is initially present for these scales. 

When a shell develops, its radius depends on the nonlinear dynamics of the chameleon 
field and is determined by Eq.(3.22) 

R^ 2 



where 



and similarly for 0oo with the background scale factor a{t) replacing R{t). 
We can estimate the size of the thin shell by rewriting: 



(5.6) 



l/(n+l) / , , . o \ l/("+l) 

0^ ^ / 3 nMpiA^ 
A r^^^^3/3fiW?(l + 5,)/ ' A V^^30:i^p^^ 



^ - I i^-(^) ...:w^'^^ J , % - At)^^^ 1 . (5.7) 



Hence we obtain 



V l/(n+l) 

^ _ 2 / MpiA \ / nMpiA^ \ 



3 , , ~ 3 



i?(t) a^(t)-i?^(t)(l + 5i)~^ . (5. 



As time evolves a{t) grows faster than R{t), which eventually stops growing and collapses 
to zero. 

We can estimate the size of the thin shell for a given model. As an example, for 
P = 0.01, n = 0.02 and A = 1.3 eV (which follows from fixing rrii = 50 Mpc^^, see 
appendix), we find approximately: 

i?2 ^ /1.6x 1020Mpc\%, , / 3,, f. :^ . x\ .. ^x 

Hence the initial thin shell disappears extremely quickly only to reappear when R is very 
small, well after virialisation. In the case of inverse power law chameleons with large 
Compton wavelength m^, we find that the thin shell is only present for a very short 
part of the cosmological evolution since matter equality. This behaviour is typical of 
inverse power law chameleon models. Of course, it would be very interesting to study 
other classes of chameleon models such as the ones inspired by f{R) models for which 
the thin shell may exist for a longer period of time. 
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5.2 f{R) models 

We can also apply our results to a class of f{R) models defined by the action 

1 



S=-^ I d^x^^fiR) (5.10) 

These models are completely equivalent to the chameleon models we have presented [24]. 
In particular, the coupling constant is uniquely fixed 

and the chameleon field is related to the curvature R via 

/'(i?) = e-2WMp, (5^12) 

while the potential is given by 

Inverse power law models can be retrieved using 

/(fi) = H_;„ + _^«(|y" (5,14) 

with -1 < p < and R:^ R:^ H^ leading to 

V^,,.,t^MkR^(-^)'^,... (,I5) 



2 p+l \ pMpi 

where Aq = '^' ° . Notice that n = — 2±i. In these models, as j3 = 1/Vq the Compton 
wave-length m^ tends to be larger than the one of the example presented in the previous 
section. Of course, a more thorough study would be necessary to determine how the thin 
shell regime could appear or disappear with more general f{R) models. 

5.3 The nonlinear regime 

We will now study the effects of the chameleon field on the evolution of spherical pertur- 
bations assuming that the inhomogeneities in the thin shell can be neglected. We first 
solve numerically Friedmann's equation for the background scale factor a{t), with time 
in units of Hubble time Hq^. We start the evolution of perturbations at a^ = 10~^. In 
the examples that follow we choose fim = 0.25, fil = 0.75, (3 = 0.01, n = 0.02 and 
A = 1.3 eV. We numerically solve Eq. (4.20) for the nonlinear evolution of the size of the 
inhomogeneity. We show in fig. 2 that the approximation mooR{a) << 1 is valid for our 
example. 

The effect of the chameleon on the evolution of the radius of the spherical perturbation 
is shown in fig. 3, where one sees that it enhances the gravitational collapse, as expected. 

19 



0.20- 



0.15 - 



0.10 



0.05 



0.0002 



0.0004 0.0006 0.0008 0.0010 



Figure 2: Plots of mooR{a) as a function of the scale factor a for Ri/ai = 100 Mpc with 
initial density contrast fixed such as to give collapse today for the usual gravity case. 
Other parameters are given in the text. 

In this example, we choose the same parameters as above and 5i = 10~'^'^^®, for which 
the collapse occurs today for the usual gravity case (/3 = 0). 

Whereas in the usual gravity case (/3 = 0) the collapse time (or acoiiapse) is independent 
of Ri, this is no longer the case in the chameleon model. Here, not only is the collapse 
time Ri dependent, but the larger the initial radius, the slower the collapse. This suggests 
that for initially large bodies, the outer shells lag behind the inner ones, never catching 
them up during collapse. Following the discussion of section 4.2, we therefore expect that 
the density of these objects will not remain uniform. Of course, the description that we 
have given is only applicable when the density is uniform in the thin shell. To gauge the 
validity of this approximation, we have tabulated in Table 1 the scale factors at collapse 
for various initial radii. We see that large enough bodies collapse with a deviation which 
is around 10 percent. Even for smaller spherical overdensities like 10 Mpc, the difference 
is only of order 25 percent. We therefore expect that our results are within the right 
ball-park. 



Ri/ai{Mpc) 


CtcoUapse 


LCDM 


1 


1 


0.528 


10 


0.769 


25 


0.823 


50 


0.857 


100 


0.886 



Table 1: Values of the scale factor at collapse for different Ri/tti 
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Figure 3: Plots of R{a) as a function of the scale factor a for Ri/ai = 1, 10 (dashed and 
dot-dashed lines) and 100 Mpc (dotted line) compared with the behaviour for R in usual 
gravity (solid line), with initial density contrast fixed such as to give collapse today for 
the usual gravity case. Other parameters are given in the text. 



Ri/ai{Mpc) 


\ia6i 


LCDM 


-3.428 


1 


-3.613 


10 


-3.497 


25 


-3.476 


50 


-3.464 


100 


-3.455 



Table 2: Values of the scale factor at collapse for different Ri/ai. 

The thin shells can be seen in fig. 4, where we plot Rg/R for Ri/ai = 1, 10,50 and 
100 Mpc, though here we have set 6i such that collapse occurs today in each case. The 
shell disappears quickly here and the subsequent evolution corresponds to the no shell 
regime. The thin shell would reappear just before collapse well after virialisation. This 
effect can also be seen in fig. 5, where we plot 7 — 1 which rapidly approaches the value 
7 — 1 = 2/3^ = 2 X 10""^ in our numerical example. The effect of modified gravity can be 
seen by in the scale dependent values of the initial density contrasts for collapse today as 
in table 2. 

Using the thin shell linear equation, let us describe the behaviour of the critical density. 
For standard gravity we find 6c{0) = 1.672, which is smaller than the usual Einstein-de 
Sitter value of Sc{0) = 1.686, since the linear growth function is suppressed in this case. 
For the chameleon case, one needs a smaller initial perturbation to reach collapse today. 
However, the linear perturbation grows faster than in standard gravity and indeed we 
find a slightly larger values of 6c corresponding to the small value of /3. For instance 
when Ri/ai = 1 Mpc, we obtain 6c = 1.677. We expect that the discrepancy would be 
stronger for models where a larger value of (3 would be compatible with a large Compton 
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Figure 4: The ratio Rs/R as a function of the scale factor a for Ri/ai 
Mpc (left to right) and the parameters given in the text. 



1,10, 50 and 100 



wavelength m'^ . 

We must emphasize that in the chameleon case, the value of the critical density is 
scale dependent, leading to a moving barrier problem for structure formation. Here this 
scale dependence is suppressed by the smallness of (3 which is numerically necessary to 
reach large values of the Compton wave length. We expect that more general models of 
the f{R) type could lead to a stronger scale dependence of the critical density contrast. 
A fully detailed phenomenological analysis is left for further work. 

6 Conclusions 

We have investigated the effects of a chameleon scalar field on the nonlinear stages of the 
gravitational collapse of a spherical overdensity region immersed in a FRW background. 
The gravitational potential is effectively changed by the presence of the new long range 
interaction. The modification depends on the profile of the scalar field. We studied the 
profile for the situation at hand, finding two very distinct classes of profiles depending 
on the size of the collapsing region. In particular, for large enough collapsing region, the 
nonlinear dynamics of the scalar field produces what is known as a thin shell effect, where 
the modification of gravity is essentially shielded. This has an infiuence on the dynamics 
of the shell as well as a potential effect on the critical density contrast for collapse, which 
in this regime of chameleon theories becomes scale dependent. 

We have also studied the spherical collapse in the case of the original chameleon model 
with an inverse power law potential and a constant coupling to Cold Dark Matter. The 
effects are larger for the smaller collapsing regions. In all cases, the initial thin shell 
disappears quickly and the collapse proceeds with no shell until a potential thin shell 
reappears after virialisation. More general models such as f{R) theories may lead to a 
stronger influence of the thin shell regime on the spherical collapse. 

We have emphasized that an inhomogeneous matter density in the thin shell will 
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Figure 5: The effect of modified gravity as measured by 7 — 1 as a function of tlie scale 
factor a and for Ri/ai = 1, 10, 50 and 100 Mpc (left to right). 

certainly result from the influence of the chameleon of the spherical dynamics. The phe- 
nomenological results in this paper have been obtained neglecting this effect. Although 
we expect that this phenomenon will be negligible for large enough initial radii, for smaller 
radii we expect that singularities such as caustics may appear. A more detailed analysis 
of these effects is left for future work. 
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Appendix: Reconstructing the Chameleon Potential 

Chameleons theory in the cosmological regime are determined by the coupling /3 and the 
time dependent mass moo(a) as a function of the scale factor as long as moo ^ H and 
the minimum of the effective potential is a dynamical attractor. Indeed, we can obtain 



V" 






ml^{a) 



/3 



2 P(q) 



(6.1: 
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Using the minimum equation we deduce 

dcj) P f> 



dt Mpi V" 

implying that 



(6.2) 



The minimum equation imphes that 

V = V, + ^ r p^-^da. (6.4) 

This defines the potential parametrically. Of course, one should check that local tests of 
gravity are then satisfied and that the long time behaviour of the model is the one of a 
A-CDM model. 

Let us give an explicit example: 

moo(a) = mi{—y (6.5) 

In the matter dominated era, we have 

V'\a) = rr^X-r - W'H^C-f (6-6) 

Assuming that rrii 3> Hi, we can neglect the second term initially. Using H{a) = 
iJi(^)3/2, we obtain 



w-*^;Si|;3T2;U-(7'^"-) <"' 



As long as 3 + 2r < 0, we can choose (pi such that 

(6.8) 



op Pi i /'^i\3+2r 



This is a power law behaviour and ^40 ^ 1 as long as f3 = 0(1). Similarly we obtain 



3/3 Vf 1 g,. 
m^M^j 6 + 2r ^ a 



v=K+i;^,rh,(-r" (6.9) 



This leads to the usual chameleon potential where we can identify 



n = -^±^ (6.10) 

3 + 2r ^ ^ 



Similarly we have 

3/3Vf 1 , 3/3p 



This is the potential we have used in the text. 
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